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I. INTRODUCTION 

Since its discovery, charmonium (bound cc states) has 
been proposed as a powerful tool to investigate many 
aspects of QCD such as the distribution of partons in 
protons and nuclei at large momentum transfer. Charm 
quarks are predominantly produced in gluon interactions 
at y/s = 200 GeV, therefore they are sensitive to the 
gluon distribution in the nucleon and its modification 
in the nucleus. In addition, the color screening of the 
cc state makes charmonium dissociation an important 
signature for the formation of a deconfined state of mat- 
ter created in A-l-A collisions [ll, 0. Such studies rely on 
a accurate understanding of charmonium production in 
p-\-p collisions which is the goal of the present work. 

The cross section of cc production is known from 
pQCD calculations to about a factor of two compared 
to PHENIX data 0, 0|- However, the hadronization step 
which forms the bound state is a nonperturbative pro- 
cess and is not well understood. A variety of schemes 
have been proposed, some of the most common being 
the Color Evaporation Model (CEM), the Color Singlet 
Model (CSM) and nonrelativistic quantum chromody- 
namics (NRQCD), also known as the Color Octet Model 
(COM). We review these models briefly here, and com- 
pare to them later in the text. 

In the CEM 0] the bound-state production mech- 
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anism is insensitive to the cc quantum numbers. The 
cc pair is produced as long as the center of mass energy 
of the pair, a/I, is greater than the mass of two charm 
quarks, but less than the mass of two open charm mesons. 
Charmonium states are then color neutralized during the 
hadronization process by soft gluon emission. The yield 
of different charmonium states is a fixed fraction J- of the 
integrated pQCD cc cross section. J- is determined from 
experiments and is universal. Hence in this model, the 
ratio between the yield of different charmonium states is 
momentum and energy independent. 



In the CSM Q the production amplitude of on-shell 
cc pairs is projected onto "^^^^Lj angular momentum 
states, and hence accounts for the J /if) and V'' as ^S'l and 
the Xc states as '^Po,i,2- The model assumes that these 
charmonium states are formed in their final color singlet 
quantum number configuration. The production density 
matrix is coupled to the wave function at the origin which 
is determined from potential models. The only empirical 
parameters entering in the entire calculation are the lep- 
tonic decay width and the charmonium mass used in the 
potential model. 



NRQCD allows for the production of both color singlet 
and color octet cc states. Color octet states emit one or 
more gluons during hadronization in order to neutralize 
their color. The production amplitude is expanded in 
powers of both the strong coupling, as, and the velocity, 
V, of the heavy quarks relative to the cc pair. The expan- 
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sion in v assumes that the heavy quark is nonrelativistic. "'^ 
As in the CSM, the production amphtudes are projected 
onto "^^^^Lj states. Since the potential model can only 
be applied to the color singlet state, a long range non- 
perturbative matrix element for each quarkonium state 
is fitted from experiments. The earliest such matrix 
parametrization W\ was tuned with J /ip and if)' cross sec- 
tions observed in CDF (^i = 1.8 TeV) P which indi- 
cated that while P-wave charmonium {xc ) has no im- 
portant color octet state contributions, S-wave charmo- 
nium (direct J /ip and ip' ) production is largely through 
color octet channels. Therefore, this model is sometimes 
simply referred to as the Color Octet Model (COM). 

Each model has its strengths and weaknesses. The 
CEM is able to reasonably describe quarkonia yields ob- 
served in many experiments, but has no predictive power 
for cc polarization. Cross sections calculated using CSM 
grossly underestimate the yields observed at PHENIX 
[10| and at CDF [§]. Recent next-to- leading order 
(NLO) [uHll and next-to-next-leading order (NNLO) 
[1Jt13] calculations for the color singlet states resulted 
in significant modifications of the predicted charmonium 
yields and polarization, but not sufficient to agree with 
the experimental results. NRQCD tuned with J/ip and 
■0' Pt spectra from CDF was able to qualitatively de- 
scribe the first PHENIX J/-0 cross section and polariza- 
tion results [l^ albeit with large experimental uncertain- 
ties, but failed to describe the J/ip and ip' polarization 
observed in CDF [l^ (see [isl for a recent review). 
Recent NRQCD calculations [20| include color singlet 
and color octet NLO short range terms along with a 
long-range matrix parametrization from experimental 
hadroproduction [2l[ and photoproduction [H, [23| of 
JJjb mesons. However, the NLO terms for the color octet 
[2J| have only small corrections compared to the leading- 
order (LO) terms and the calculations still disagree with 
the J/-0 polarization measured by CDF. 

One of the complications in the J/ ip total cross section 
and polarization calculations (observables where exper- 
imental tests are readily available) is the contribution 
from the decays of excited charmonium states, primarily 
Ip' , Xci and Xc2- In addition, the J/ip suppression ob- 
served in heavy ion collisions cannot be completely un- 
derstood without a knowledge of the feed-down fraction 
of excited charmonium state decays to the J/ip. This is 
particularly true under the assumption that the suppres- 
sion is due to the disassociation of charmonium in the a 
high temperature quark-gluon plasma, since lattice calcu- 
lations [2] indicate that the melting points of the Xc and 
ip' states are lower than that of the J/tp . In this work 
the feed-down fractions to the J/ tp from excited charmo- 
nium states is measured since they can be determined 
more precisely than production cross sections as many 
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of the systematic uncertainties cancel when making cross 
section ratios. 

The PHENIX Experiment at RHIC can measure 
quarkonia dilepton decays over a broad pr and rapidity 
range and can detect photons from Xc radiative decays 
using electromagnetic calorimeters at midrapidity. This 
article reports on the feed-down fraction of the J/tp yield 
which comes from ■(/;' and Xc decays at midrapidity in 
p+p collisions at y/s = 200 GeV. For these we used the 
i/y to J/ip yield ratio in the dielectron channel and the 
full reconstruction of the p+p Xc J/ip + j ^ 
e+e~-|-7 decay. A new J/ip differential cross section mea- 
surement at midrapidity and forward rapidity using the 
increased luminosity obtained in the 2006 and 2008 runs, 
is also presented. These provide more accurate measure- 
ments than previously published in jTol |. particularly for 
the J/ip differential cross section at high pr ■ The re- 
sults obtained in these analyses also provide a baseline for 
the study of J/ip suppression in d-|-Au [l^ and Au-|-Au 
[ii, [23 colhsions at PHENIX. 

Systematic uncertainties throughout this article are 
classified according to whether or not there are point- 
to-point correlations between the uncertainties. Type A 
systematic uncertainties are point-to-point uncorrelated, 
similar to statistical uncertainties, since the points fluc- 
tuate randomly with respect to each other. Type B sys- 
tematic uncertainties are point-to-point correlated. The 
points fluctuate coherently with respect to each other. 
That is, it accounts for the uncertainty in the nth-order 
derivative of the measured spectrum, in most cases the 
slope. Global, or type C, systematics are those where all 
points fluctuate in the same direction and by the same 
fractional amount. 

The remainder of the article is arranged as follows. 
An introduction to the PHENIX detector, a description 
of the data sample, and a description of the lepton iden- 
tification method is described in Sec. [Hi The analy- 
sis is described in three sections: midrapidity J/ip and 
ip' dielectron cross section measurement in the PHENIX 
central arms (Sec. IIIip : direct Xc feed-down measure- 
ment in the central arms (Sec. IIVP : and forward rapidity 
J/ip cross section measurement in the muon arms (Sec. 
\V\ . The results are compared to measurements from 
other experiments and to current theoretical calculations 
in Sec. IID 



II. EXPERIMENTAL APPARATUS AND THE 
DATA SET 

The PHENIX detector [28| is composed of four 
arms. Two central arms measure electrons, photons and 
hadrons over \r]\ < 0.35 with each azimuthally covering 
A(p — Tr/2. Two forward muon arms measure muons over 
the range —2.2 < rj < —1.2 arm and 1.2 < ry < 2.4 with 
full azimuthal coverage. 

Charged particle tracks in the central arms are formed 
using the Drift Chamber (DCH), the Pad Chamber 
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(PC) and the collision point. Electron candidates 
required at least one fired phototube within an an- 
nulus 3.4 < Rring[cm] < 8.4 centered on the projected 
track position on the Ring Imaging Cerenkov detec- 
tor (RICH).^ In addition, the electron candidate is re- 
quired to be associated with an energy cluster in the 
Electromagnetic Calorimeter (EMCal) that falls within 
^cfposition of the projected track position, and within 
4^'^ energy of the cxpected encrgy/momentum ratio where 
the ct's characterize the position and energy resolution 
of the EMCal. The relatively loose association require- 
ment still provides excellent hadron rejection due to the 
very small particle multiplicity in p+p collisions. Based 
on the pt range of decay electrons from J/^p observed in 
real data and simulations, a minimum px of 500 MeV/c 
was also required for each electron candidate. 

Each of the forward muon arms [2^ comprises a hadron 
absorber, three stations of cathode strip chambers for 
particle tracking (MuTr), and a Muon Identifier detector 
(MuID). The hadron absorber is composed of a 20 cm 
thick copper nosecone and 60 cm of iron which is part of 
the magnet. The MuTr is installed in an eight-sided con- 
ical magnet. The MuID is composed of five steel hadron 
absorbers interleaved with six panels of vertical and hori- 
zontal larocci tubes. A single muon needs a longitudinal 
momentum of 2 GeV/c to reach the most downstream 
MuID plane. Tracks reconstructed in the MuTr are iden- 
tified as muons if they match a "road" formed by hits 
in the MuID, within 2.5cr of angular resolution. At least 
one tube in the last MuID plane should have fired. Addi- 
tional cuts include a < 23 for the reconstructed track, 
a < 9 for the track projection to the collision vertex 
and a polar angle cut of 14° < 0^ < 33° for the north 
arm and 147° < 0^ < 166° for the south arm to avoid ac- 
ceptance inconsistencies between the detector simulation 
and real data near the edges of the muon arms. 

Beam interactions were selected with a minimum-bias 
(MB) trigger that requires at least one hit per beam 
crossing in each of the two beam-beam counters (BBC) 
placed at 3.0 < |?7| < 3.9. Studies using Vernier 
scans (also called van der Meer scans) [10] conclude 
that this MB trigger accepts a cross section of ctbbc = 
23.0 ± 2.2 mb. This cross section represents 55 ±5% 
of the cr™*^'' = 42 ± 3 mb p-\-p inelastic cross section at 
= 200 GeV. 

Dedicated triggers were used to select events with at 
least one electron or two muon candidates. An EM- 
Cal RICH Trigger (ERT) required a minimum energy 
in any 2x2 group of EMCal towers, corresponding to 
A77 X = 0.02 X 0.02 rad., plus associated hits in the 
RICH. The minimum EMCal energy requirement was 400 
MeV for the first half of the Run and 600 MeV for the sec- 
ond half. The data used in this analysis were taken with 
the ERT in coincidence with the MB trigger. Events were 



also triggered when there were two muon candidates in 
one of the MuID arms. The trigger logic for a muon can- 
didate required a "road" of fired larocci tubes in at least 
four planes, including in the most downstream plane rel- 
ative to the collision point. The event sample used in the 
dimuon analysis required a MuID trigger in coincidence 
with the MB trigger. 

There are events which produce a J but do not fire 
the MB trigger. The fraction, Sinei.^oi such events is es- 
timated by measuring the number of high pT 7r° — >■ 77 
decays which satisfy the minimum energy condition of 
the ERT and which do not satisfy the MB trigger. It 
was found that Einei. = (79 ± 2)%. The correction due to 
this factor is included in all cross section calculations for 
measurements requiring the MB trigger. No dependence 
of Sinei. on the Pt of the measured 7r° decays was found 
over the range 0-10 GeV/c [H. 

The collision point along the beam direction is deter- 
mined with a resolution of 1.5 cm by using the difference 
between the fastest time signals measured in the north 
and south BBC detectors. The collision point was re- 
quired to be within ±30 cm of the nominal center of 
the detector. In the dielectron analysis, runs in which 
electron yields were more than three standard deviations 
away from the average in at least one of the eight EM- 
Cal azimuthal sectors, were discarded. For the dimuon 
analysis, runs where the muon arm spectrometers were 
not fully operational were rejected. 

The 2006 data sample used in the dielectron analysis 
corresponded to Npp = 143 billion minimum bias events, 
or an integrated luminosity of J Cdt — Npp/aQQc = 
(6.2 ± 0.6) pb-^ The 2006 and 2008 data samples used 
for the muon analysis, corresponded to 215 billion mini- 
mum bias events, or a luminosity of (9.3 ± 0.9) pb^^. 



III. J/il> AND tp' ANALYSIS IN THE 
MIDRAPIDITY REGION 

The procedure for analyzing the J/ip and ip' — > dielec- 
tron signal in the central arm detectors is detailed in this 
section. The overall procedure to select dielectrons and 
extract the charmonium signal and determine combina- 
torial and correlated backgrounds is explained in IIII Al 
Studies of the central arm detector response to charmo- 
nium dielectron decays is the subject of the section UlIB I 
The final pt and rapidity dependence of the cross sections 
is calculated in Section IIII CI together with a summary 
of all systematic uncertainties mentioned throughout the 
text. Finally the ip' / {J / -ip ) dielectron yield ratio is cal- 
culated in Section IIIIDI 



A. Di-electron decays of J/xp and %p' mesons in the 
midrapidity region. 



^ Corresponding to A<?!> = 8mrad and AZ = 3cm 



The invariant mass was calculated for all electron pairs 
in which one electron of the pair geometrically matched 
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ee invariant mass [GeV/c ] 

FIG. 1: (Color online) Invariant mass distribution of unlike- 
sign (closed circles) and like-sign (open boxes) dielectrons in 
the J/i/) and tp' region without (a) and with (b) a minimum 
Pt requirement for the dielectron pair. Dash-dotted (dashed) 
lines represent the mass range used to count J /ip {ip' ) decays. 



the position of a fired ERT segment. This requirement 
was necessary given that we used simulated J/ip and 
tp' decays to estimate the ERT efficiency. Di-electron 
contributions to J/ip and iJj' decays are clearly identified 
as peaks in this invariant mass distribution (Fig. [ij. 
The primary sources of physically correlated unlike-sign 
pairs (e^e~) are quarkonia decays, open cc and 66 pairs, 
Drell-Yan, and jets. Uncorrelated unlike-sign pairs are 
from combinatorial background. The primary sources 
of like-sign pairs (e^^e^ -|- e~e~) are combinatorial back- 
ground, and electrons from particle decays occurring in 
the same jet (mostly 7r° Dalitz decays). The like-sign 
pair mass distribution normalized by the geometric mean 
of the number of e"'"e"'" and e~e~ pairs was statisti- 
cally subtracted from the unlike-sign mass distribution. 
The primary effect of this subtraction was to account 
for combinatorial background, however it also accounted 
for much of the jet background. There are 2882 unlike- 
sign and 203 like-sign dielectrons in the J/^ mass range 
(2.7< Mee[GeV/c^ ] <3.4), giving a correlated signal of 
2,679 ± 56 counts and a signal/background of 13. In the 
-0' mass region (3.5< Mee[GeV/c^ ] <3.9) there were 137 
unlike-sign and 51 like-sign electron pairs corresponding 
to a signal of 86 ± 14 counts and signal/background of 
1.7. 

The jet contribution in the charmonium mass region is 
three orders of magnitude smaller than from the J/ip and 
ip' with a steeply falling mass spectrum[32:] and will be 
ignored here; in any case it is largely removed by the 



like-sign subtraction. The Drell-Yan contribution was 
estimated using next-to- leading-order calculations [33 |. 
Taking into account the detector acceptance, the frac- 
tion of the dielectron signal which comes from Drell-Yan 
processes is 0.23 ± 0.03 % in the J/ip mass region and 
3.37 ± 0.40 % in the ip' region. The heavy quark con- 
tribution is the major background to the correlated di- 
electron spectrum. In fact, they represent a significant 
fraction of the correlated dielectrons in the tp' mass re- 
gion. They will be estimated by two models as described 
in the next several paragraphs. 

In order to understand the dielectron spectrum, a sim- 
ulation was done for the three primary contributions to 
the mass spectrum: the J/ip and ?/;', heavy quark pairs, 
and Drell-Yan. The first step was to generate the initial 
correlated electron pair spectrum. The J/ip and ip' were 
generated by weighting their distributions in order to ob- 
tain the same px spectrum as seen in real data. The 
J/ip radiative decay (J/V' e~^e~ also called in- 

ternal radiation, was introduced using the mass distri- 
bution estimated from QED calculations [sij. Drell-Yan 
pairs were generated according to the mass distribution 
obtained from NLO calculations. In order to make a 
conservative estimate and determine whether the result 
is model independent, the cc and bb mass distributions 
were obtained using two different methods: 

1. A dielectron generator: The semi-leptonic 
heavy flavor yield measured in [sil . Issj was split into the 
cc {dcTcc/ dpx) and bb [da^i/ dpT) distributions according 
to the c/b ratio from fixed-order plus next-to- leading- 
log (FONLL) calculations [1] which agree with PHENIX 
measurements of separated cc and 66 production [36l [3^ . 
These cc and 66 yields were used as input for an electron 
Monte Carlo generator with uniform rapidity distribu- 
tion {\y\ < 0.5) and the measured vertex distribution. An 
electron and positron from the decay of a heavy quark 
pair were generated for each event. In this method the 
heavy quarks are assumed to have no angular correlation. 

2. PYTHIA: Hard scattering collisions were simulated 
using the PYTHIa[3^ generator. Leading order pair cre- 
ation sub-processes and next-to-leading-order flavor cre- 
ation and gluon splitting sub-processes are all included 
in the heavy quark generation 39]. These sub-processes 
have different opening angles for the heavy quark pair. 
The simulation used the CTEQ6M parton distribu- 
tion functions (PDF), a Gaussian fcy = distribution of 
width 1.5 GeV/c , a charm quark mass of 1.5 GeV and 
bottom quark mass of 4.8 GeV. Variations of the fcy dis- 
tribution and masses of the heavy quarks were included 
in the systematic uncertainties. The pt dependence of 
electrons from cc and 66 given by the simulation agrees 
with the PHENIX measurement of single electrons from 
heavy flavor decavfsij. 

The generated electron pairs from all sources were then 
used as input to a GEANT-3 [4l| based detector Monte- 
Carlo which included effects such as Bremsstrahlung ra- 
diation of electrons when crossing detector material and 
air (external radiation). Simulated events were then re- 
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constructed and analyzed using the same criteria as were 
used for real data and reported in Sec|lT] and Sec lIII Al 
More details will be given later in Sec lIII B| including 
methods of estimating systematic errors. 

The resulting simulated distributions were then 
fit to the 2-dimensional mass vs pt distribution 
of the measured dielectron signal in the mass 
range 2.0 < M^+^-[GeV/c^ ] < 8.0. The Jit parame- 
ters included the normalization of cc , bb , J/ip and 
^p' contributions, the fraction of the internal radiation, 
and a mass resolution correction for the simulated res- 
onance peaks. The normalization of the Drell Yan was 
fixed according to expectations from the NLO calcula- 
tions. 

Fig. [5] shows the results of the fit for the dielectron 
mass (a) and px (b,c). The heavy flavor contribution to 
the continuum obtained from the fit using the dielectron 
generator and pythia is shown in Fig. [2]-d. When us- 
ing the PYTHIA simulation, the presence of back-to-back 
correlated cc and bb pairs produced more high mass pairs 
per cc which then forced a smaller contribution from bb. 
As can be seen from the figure, the fits performed us- 
ing the two generators give very different normalizations 
for the open charm and the open bottom contributions. 
However, the two methods give very similar contributions 
for the sum which is well constrained by data. Thus 
the lack of the knowledge of the angular correlation in 
heavy flavor production does not affect the estimate of 
the total continuum contribution from open heavy fla- 
vor in the J/ip and tp' mass regions. The measurement 
of the cc and bb cross sections is not in the scope of this 
paper; a more detailed study can be found in j37ll42l. |43| . 
Type A fit parameter uncertainties and the type B un- 
certainty obtained from the difference in results obtained 
using the two generators for the total heavy fiavor contri- 
bution, are summed in quadrature and shown as bands 
in Fig. [2] Values for the fraction of the charmonium sig- 
nal (/^) shown in Figs. [5]-b and[5]-c are used later in the 
yield calculation. 

The fitted external and internal radiation contribu- 
tions indicate that the fraction of radiative decays of 
the J/ip , where the undetected photon has energy 
larger than 100 MeV, is (9 ± 5)%. This is consis- 
tent with QED calculations which indicate that 10.4% 
of the dielectron decays from the J/ip come from such 
radiative decays and a measurement of fully recon- 
structed J/ip e+e^7 performed by E760 44] which 
gives 14.7 ± 2.2 %. The J/i/j mass peak around 
3.096 GeV/c^ has a Gaussian width from the fit of 
53 ± 4 MeV after including a mass resolution in the 
MC of (SM/M) of (1.71 ± 0.13)%. Because of the ra- 
diative tails, the mass range (2.7< Mee[GeV/c^ ] <3.4) 

contains Smass = 93.8 ± 0.9% of the Jf-tp decays and 
the mass region (3.5< Afee[GeV/c^ ] <3.9) contains 



^mass = 86 ± 2% of the 'ip' decays, corrections included 
in the yield calculations. The foreground yield as well as 
the statistical uncertainties used in the cross section cal- 
culations were obtained assuming that both foreground 



and background distributions are independent and follow 
Poisson statistics. The total foreground was then multi- 
plied by the factors obtained from fits in the previous 
section to obtain the J/ip and ip' yields. In each bin of 
Pt (or y) the foreground signal (/i/) was obtained from 
the unlike-sign counts (fg) in the distribution and the 
background (fib) was obtained from the like-sign counts 
(bg) (Fig. [T]top). The joint probability distribution for 
the net number of counts s — fif — fif, 



We expand the term + ■ 



(1) 



s \ 



fg fa 



l + — 



' k=0 



ifg-ky.kl VMb 



(2) 



^^''^'^ = lg\{fg-k)\ fcT 



fc=0 



Assuming no negative signal, the expression is summed 
over from to 00 using the normalization of the 
Gamma distribution 



dxx^ e 



i„-bx _ P - 1 



bP 



(3) 



and b = 2, p — 1 = bg + fg ~ k. We obtain finally. 



fa 



{bg + fg-k)n (I 



k=0 



bg\{fg-k)\) 2 V2 



bg+fg-k gfcg-s 



fc! 



(4) 



The number of charmonium decays for each bin, and 
the corresponding statistical uncertainty, were obtained 
using Q given the fraction (/^) of charmonium in the 
sample found previously. 



(5) 



B. Di-electron acceptance and efficiency studies 

The detector response to J/ip and ip' dielectron decays 
was studied using the GEANT-3 based Monte Carlo sim- 
ulation. Malfunctioning detector channels were removed 
from the detector simulation and from the real data anal- 
ysis. The geometric acceptance of the detector Monte 
Carlo was compared to that for real data using simu- 
lated 7r° decays. The majority of the electrons found in 
real data come from tt^ Dalitz decays and photons which 
convert to electrons in the detector structure. The simu- 
lated electrons from tt^ decays were weighted in order to 
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FIG. 2: (Color online) Correlated dielectron mass (a) and pr distributions in the J/tp (b) and ?/>' (c) mass regions. Signal com- 
ponents were estimated by fitting open heavy flavor, Drell Yan (normalization fixed by NLO calculations), J/tp and t/i' decays 
after detector simulation. The cc and bb components were generated using pythia [3g| and a heavy flavor based dielectron 
generator described in the text. Bands correspond to the type A fitting uncertainties and the type B systematic uncertainty 
obtained when using the two different open heavy quark generators. Panel (d) shows the result after the fit for cc, bb and total 
open heavy flavor components from each generator. 
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FIG. 3: (Color online) Simulated (solid line) and real data 
(points) single electron distributions in the <j) coordinate of 
the drift chamber. Error bars correspond to statistical uncer- 
tainties. 



match the coUision vertex and pt distributions observed 
in the data. Fig. [3]shows the simulated and real electron 
track distribution as a function of the azimuthal angle, 
(f>, measured at the DCH radius. The ratio between real 
and simulated track distributions {facc{4>DCH, zjjch)) is 
used later to estimate the systematic uncertainty of the 
J/ip acceptance. 

The electron identification efficiency was estimated us- 
ing 7 — !■ e+e^ conversions coming primarily from the 
beam pipe. These dielectrons, which do not originate 
from the event vertex, have a nonzero invariant mass 
and can be identified since their invariant mass exhibits a 
peak in the region below 30 MeV/c^. Assuming all tracks 
in the peak above the combinatorial background are elec- 
trons, the electron identification efficiency was obtained 
from the fraction of dielectron conversions which survive 
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FIG. 4: Single electron identification efficiency estimated us- 
ing photon conversions from real data (points) and the elec- 
tron simulation (shaded area). 
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FIG. 5: Transverse momentum dependence of the detector 
performance and correction factors for dielectron decays of 
J/t/j mesons in |y| < 0.5. Shaded bands are the uncertainties 
of the estimates as described in the text. 



the identification criteria applied to both electron and 
positron compared to the number of dielectron conver- 
sions obtained after requiring identification for only one 
electron or positron. The same procedure was repeated 
in the simulation. Fig. |4] shows the electron identifi- 
cation efficiency as a function of the pr of the electron 
in question. The difference in efficiency between simu- 
lation and data for electrons with pT > 0.5 GeV/c was 
no larger than 0.8%, which translated to an overall type 
B uncertainty in the dielectron yield of 1.1% due to our 
understanding of the electron identification efficiency. 
Simulated J/ip dielectron decays were generated with 



uniform pT and rapidity (\y\ < 0.5) and the mea- 
sured vertex distribution. The fraction of the generated 
J ftp decays that were fully reconstructed corresponds to 
the acceptance x electron identification efficiency of the 
detector {Ax Eeio) for J/'ip dielectron decays with rapid- 
ity \y\ < 0.5 (Fig. [5]-a). When each simulated electron de- 
cay was weighted according to facc{4>DCH,ZDCH) given 
previously, the number of reconstructed J /tp decays was 
modified by 7.5%. This is essentially the variation in 
our acceptance calculation, when calculated using a data 
driven method as compared to simulation. We considered 
this deviation as a type B systematic uncertainty. The 
A X SeiD for simulated ip' dielectron decays in the same 
rapidity range was larger than that from the J/ip by be- 
tween 5-20% because of its larger mass. The maximum 
difference occurs at px ~2.5 GeV/c . 

The detector acceptance for charmonium also depends 
on the orientation of its electron decay with respect to 
the momentum direction of the parent particle, an out- 
come of charmonium polarization. The correction factor 
from polarization (Spoi) was evaluated using a measure- 
ment of J/ip polarization [l8| in p+p collisions interpo- 
lated to the relevant transverse momentum. The uncer- 
tainty in Epoi due to the uncertainty in the polarization 
was assigned as a type B systematic uncertainty. In the 
Pt region where there is no polarization measurement 
[pt > 5 GeV/c for J/4> and all pr for ijj' ) the one stan- 
dard deviation uncertainty was calculated assuming the 
J/V' polarization in this region can be anything between 
-1 and 1. Fig. [5]-b shows the pr dependence of Spoi- 

The trigger (ERT) performance was studied using sin- 
gle electrons. We used a MB data sample to measure 
the Pt dependent fraction of electron candidates that 
fired the ERT in each of the EMCal sectors. These 
fractions were then used in simulation to estimate the 
J/ip efficiency of the ERT trigger (eert)- This process 
was repeated for each change in the ERT operational con- 
ditions, such as a change in the energy threshold, or a sig- 
nificant modification in the number of EMCal or RICH 
sectors included in the ERT trigger. Fig. [5]-c shows the 
Pt dependence of Ser.t, weighted by the luminosity ac- 
cumulated in each ERT period. When the single electron 
ERT efficiency of each EMCal sector was varied within its 
statistical uncertainty, a one standard deviation change 
of 4.5% in Eert was observed. This deviation is shown 
in Fig. [5]-c as the shaded band and is assigned as a type 
B systematic uncertainty for the J/ip and ip' yields. No 
significant change in Eert was observed if one used the 
-0' in the simulations. 

A final correction {bin shift) was made for the domi- 
nance of the yield in the lower end of each pT bin (Fig. 
[S]-d). In addition, a correction of up to 2% was made 
to account for bin-by-bin smearing effects due to finite 
momentimr resolution(6insmear) • 
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C. Cross section results 

The J /'ip and ip' dilepton differential cross section for 
each pt bin is calculated by 



Bt d'<7. 



2TrpT dydpT 



B. 



1 d^N/dydpT 

2tTPt Einel. J Cdt 

N,i, 



(6) 



dprdy AyApxAe 

is the branching ratio of 
states into dileptons 



the 
and 



where 

charmonium 

smear ■ 

All systematic uncertainties described in the previous 
sections are listed and classified in TablelH The quadratic 
sum of the correlated systematic uncertainties (type B) is 
between 10% and 13% of the measured J ftp yield and be- 
tween 12% and 22% of the measured ip' yield, depending 
on pt ■ 



TABLE I: List of the systematic uncertainties relative to 
the J ftp and ip' dielectron yields. Ranges indicate pr depen- 
dence. 



description 


contribution 


type 


fraction of J/i/' in the mass cut 
fraction of ip' in the mass cut 


0.4% 


A 


3-13% 


A 


acceptance 


7.5% 


B 


elD efficiency 


1.1% 


B 


mass cut efficiency for J/ip 
mass cut efficiency for 


1.0% 


B 


2.0% 


B 


heavy flavor MC used in fit for J/tp 
heavy flavor MC used in fit for ^' 


0.5-1.1% 


B 


4.8-10% 


B 


up-in-down bin correction 


3% 


B 


momentum smear effect 


1.5% 


B 


Pt , y and vertex input in xp MC 


2.0% 


B 


J/tp polarization bias in acceptance 
ip' polarization bias in acceptance 


0-10% 


B 


4-17% 


B 


ERT efficiency 


4.5% 


B 


luminosity 


10% 


C 



The Pt dependencies of the measured J/^p and 
ip' yields are shown in Fig. EJtop) and Tables I VIII 
IVIIII The bars correspond to the quadratic sum of all 
type A and statistical uncertainties. Boxes represent the 
quadratic sum of the type B uncertainties. There is a 
global uncertainty (type C) of 10%. 

The PT integrated J/^ cross section was calculated for 
three rapidity ranges using 



" dy ^ "dpTdy 



Apj 



(7) 



where Bud^a / dpxdy is obtained from (|6]) using Aseio, 
£ert E^nd Epoi recalculated for each of the three rapidity 
bins. The results are listed in Table HXl and shown in Fig. 

m 
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FIG. 6: Transverse momentum dependence of J/^p and 
tp' yields in \y\ < 0.35 (top), tp' /{J/ip ) ratio together with 
ratios obtained in other experiments (bottom). Error bars 
reflect statistical and type A uncertainties while boxes reflect 
the quadratic sum of type B uncertainties). 



D. ip'/(J/^p ) yield ratio and fraction of J/^p yield 
coming from ^p' decays. 

The decay of ^' to J/tp cannot be measured in the cur- 
rent detector configuration. However, we can calculate 

the fraction of J/ip coming from ip' decays ^^^^/^^'^ using 

the ratio between the ip' and J/ip cross sections and the 

iP' branching ratio to J/ip (B^/^ = (58.7 ± 0.8)% 



(8) 



We start from the ratio between the ip' and the 
J/tp dielectron counts R^/^- Its joint probability distri- 
bution is calculated from the expected Poisson probabil- 
ity distributions (Eq. |4]) P^'{s^') and Pj/^{sj/^) for the 
dielectron counts in the ip' and J/ip mass ranges respec- 
tively, and the corresponding values and fj/^' which 
account for the fraction of ip' and J/ip contributions in 
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the chosen dielectron mass ranges: 



Pi? 



4' 



(9) 



The ijj' /{J/tfj ) dielectron cross section ratio is thus de- 
termined as follows where the different correction factors 
for ^p' and J/ip must be taken into account. 



^ R 



_ , , j/i> j/i> 



(10) 



Type A uncertainties are propagated for and fj/^' 
while common relative type B uncertainties that are cor- 
related for J/ip and ip' cancel. The remaining uncer- 
tainty in the ratio comes from the quadratic difference 
between type B uncertainties which are different for the 
J/ip and ip' . The ip' /{J/ip ) dielectron cross section ratio 
is shown in the bottom panel of Fig. [B] The numbers are 
listed in Table H 

Using the branching ratios, B^+^_ = (0.765 ± 0.017)% 

and B^l'l!_ = (5.94 ± 0.06)% gll in ® gives 



F^/'^ = (9.6 ± 2.4)%. 



IV. RADIATIVE DECAY OF xc 



(11) 



The decay channel Xc — ^ J/fp +7— >-e+e^ +7 is fully 
reconstructed in the central arms and is used to directly 
measure the feed-down fraction of Xc decays in the inclu- 
sive J/ip yield {F^i^)- This measurement is particularly 
challenging since the photon is typically of very low en- 
ergy. The data sample used in this measurement and the 
7 identification procedure is described in Section IIV Al 
The detector performance for the measurement of photon 
decays of the Xc is discussed in Section HVBI The com- 
position of all combinatorial and correlated backgrounds 
for the Xc signal in the e"'"e~ 7 mass distribution is de- 
tailed in Section IIV CI Section IIVDI presents the final 
feed-down fraction calculation and a summary of all un- 
certainties. 



A. Selection of Xc ^ J /''P + 7 decays 

The analysis of the radiative decay of the Xc requires 
the identification of photons with energy {E^) as low as 
300 MeV, the lower limit of the energy we allow in this 
analysis. Photons were identified as energy clusters in 
the EMCal whose profile is consistent with an electro- 
magnetic shower. This profile is based on the response of 
the EMCal to electron beam tests performed before the 
EMCal installation . Energy clusters that were closer 
than four standard deviations (of the energy cluster po- 
sition resolution) to reconstructed charged tracks were 
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FIG. 7: (Color online)Study of the tt" detection performance 
using both the measured 7 energy in real data (boxes) and 
detector MC (shaded band), (a) 7 identification efficiency, (b) 
tt'' mass peak position, (c) tt" mass resolution. The vertical 
dashed line represents the minimum 7 energy required in the 
Xc analysis. Uncertainties are from the vr" fit parameters in 
simulated and real data. 



rejected, in order to remove electron and misidentified 
hadron contributions. Electrons from photon conversions 
in detector material which were not reconstructed by the 
tracking system were removed by requiring energy clus- 
ters to be further than four standard deviations from hits 
in the Pad Chamber (PC) located in front of the EMCal. 

The invariant mass of e^e'^ + 7 is formed us- 
ing e+e^ pairs in a tight J /%p mass region of 2.9 < 
A4+g_ [GeV/c^ ] < 3.3, avoiding the region where pho- 
tons produced by Bremsstrahlung radiation can become 
an additional background in the 300 MeV energy region. 
The sample contains Nj/^ =2456 ± 51 e^e~ pairs from 
J decays, after removing combinatorial and correlated 
background as done previously. The e"*"e~7 mass dis- 
tribution is plotted in Fig. [T^] (top), where we require 
E^ >300 MeV. The mass of e+e^ 7 minus the mass of 
the measured e^e~ pair is plotted in order to cancel the 
effect of the mass resolution in the e+e~ pair. The re- 
maining resolution in the subtracted mass distribution is 
from the energy resolution of the measured photon. 



B. Detector performance for Xc radiative decay 



The resolution of the mass distribution Mg+g-^ — 
Mg+g- is dominated by the photon energy resolution of 
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the EMCal. Most photons from Xc decays have energy 
close to the lower limit of the EMCal sensitivity. The 
behavior of the calorimeter was studied by using a clean 
sample of tt*^ — >■ 77 decays in real data and in the simula- 
tions. Pairs of clusters were formed where the invariant 
mass of the pair was required to be consistent with a 
tt" . Only one of the clusters was required to pass electro- 
magnetic shower requirements. The photon identification 
efficiency was obtained assuming the other cluster of the 
pair was a photon. This was done on a statistical ba- 
sis by subtracting a mixed event background to account 
for the small contamination from random clusters under 
the 7r° peak. Fig. [7]-a shows the energy dependence of 
the photon identification efficiency (e^/d) obtained using 
real and simulated tt'^s. The simulation gives an efficiency 
2.3% larger than that found in real data. This difference 
was assigned as a type B systematic uncertainty in Sjid- 
The central value of the tt'^ mass peak decreases slightly 
as the photon energy approaches the lower limit of the 
calorimeter sensitivity. This behavior is caused by zero 
suppression during data acquisition and the energy clus- 
ter recognition algorithm. These effects are correctly re- 
produced in simulation as can be seen in Fig. [7]-b. The 7 
energy resolution {5Ej/Ej) was uniformly degraded by 
4.7% in the simulation in order to match the mass res- 
olution (SM/M) of the 7r° peaks observed in real data 
(Fig. He). 
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FIG. 8; (Color online) Mass (a) and 7 energy (b) distribu- 
tions of e'*'e~ 7 decays from xa and Xc2 decays obtained from 
PYTHIA -I- detector simulation. The dashed line in (b) repre- 
sents the 7 energy cut applied in this analysis. 

Xci and Xc2 states were generated using gluon-|-gluon 
scattering in PYTHIA with the CTEQ6M PDF, requiring 
that the J/^ be in the rapidity range \y\ < 0.5. The Xco 
is not considered in the simulation because of its small 
branching ratio to J/V- of (1.14 ± 0.08)% Fig. [8] 

shows the mass and 7 energy distribution of e"'"e~ 7 de- 
cays of simulated Xc- The conditional acceptance of 7 
from Xc is plotted as a function of the J/tp momentum 
in Fig. |9l The detector geometric acceptance of the 
Xc can be affected by its polarization and the polariza- 
tion of the decay J/tp . There is no measurement of the 
Xc polarization. Simulation studies found the overall ac- 



ceptance is modified by at most 5.6% if the Xc is totally 
transversely polarized. This possible modification was in- 
cluded in the acceptance type B systematic uncertainty. 




^ P^(J/V) [GeV/c] 



FIG. 9: Conditional acceptance and efficiency of the Xc decay 
7 as a function of the J/ip transverse momentum. The height 
of the boxes corresponds to the type B systematic uncertainty 
due to the lack of knowledge of the polarization and the pho- 
ton identification efficiency. 



C. Composition of the 



+ 7 sample. 



In addition to the Xc signal, the observed e+e^ -f- 7 
sample is composed of combinatorial background, mostly 
coming from uncorrelated tt" decays present in events 
where a J/ip is detected, and by photonic sources corre- 
lated to the J/ip which will be discussed later. 

The combinatorial background from random 
e+e^ pairs (i.e. the combinatorial background to 
the J/ip in the Xc decay) is well described by the sum 
of e+e~'"7 and e~e~7 mass distributions. This sum 
was normalized by the geometrical average of the two 
components, and subtracted from the e'^e~ 7 mass 
spectrum. The mass distribution of random (e+e^ ) -I- 7 
combinations (i.e. essentially random J/ip -I-7 pairs) 
was obtained using the invariant mass distribution of 
e+e~ pairs from one event and photons another. In 
order to obtain the combinatorial background as realis- 
tically as possible, events used to form the e+e" and 7 
combination were required to have event vertices within 
3 cm (2(T of the vertex position resolution) of each other. 

The sources of correlated background include internal 
and external (Bremsstrahlung) radiative decays of J/ip , 
i.e. J/ip — > e+e^ 7, tt'^s produced in jets containing 
J/ip , V' ~^ J/i' + neutral mesons, i?" — J/4> + ^ 
where X or its decays includes a 7. Another possibil- 
ity is that a J/tp could be produced together with a high 
energy photon[47|. Recent studies also suggest an impor- 
tant contribution from gg — >■ J/ip +'^ + gg in NNLO cal- 
culations at y/s =14 TeV [l6|. No estimate was made for 
y/s — 200 GeV at the time of this writing. These sources 
will be considered in the next few paragraphs. 

Photons produced by Bremsstrahlung radiation in the 
detector structure are very close to their associated 
electron and are rejected by the criteria that removes 
electrons in the 7 identification. The minimum di- 



13 



electron mass cut of 2.9 GeV/c^ also removes radiative 
J /ip decays with > 200 MeV, i.e. those in the energy 
range of the photons used in this analysis. 

Collisions containing primary J/^ mesons produced 
by gluon+gluon scattering (the dominant source) were 
simulated using pythia in order to understand the elec- 
tron radiation and jet contributions. Only the e'^e~ and 
the radiative e~^e~ 7 decay channels were allowed. All fi- 
nal state particles with momentum larger than 100 MeV 
and I77I < 0.5 were reconstructed. J/tp and 7 identifi- 
cation criteria were the same as used in the analysis of 
real data. The e"'"e~7 distribution obtained from this 
simulation is completely accounted for by combinatorial 
background from mixed events (Fig. rTOf a)). leaving lit- 
tle room for contributions from possible jets containing 
J/ijj , radiative decays or electron radiation when cross- 
ing the detector support. 

Using the data, a check was done for possible miss- 
ing correlated radiation backgrounds that might have 
been missing in the simulation. The invariant e'^e~j 
mass distribution was formed in which we required 
Mg+e- [GcV/c^ ] <2.9. The Xc contribution is small in 
this region and the correlated signal should be mainly 
from other sources, e.g. J/ip internal and external radi- 
ation. The data unlike the simulation shows a correlated 
background after combinatorial background subtraction 
(Fig. [TOl b)). The line shape of this mass distribution can 
be described by a Gaussian distribution, Landau distri- 
bution, or a simulated ip' ^ J + 7 shape. Its source 
could be the gg ^ J/ip +j + gg process mentioned pre- 
viously, but we simply take this as a background which 
must be included in the fit to the Xc mass distribution. 
The position of the peak is set by the minimum pho- 
ton energy cut of 300 MeV, while the width is set by 
energy spectrum of the source and more importantly by 
the smearing effect caused by the fact that the spectrum 
is a difference of two invariant mass calculations. These 
results will be used later when fitting the e"'"e~ 7 invari- 
ant mass distribution. 

In section HITD] we reported that (9.6 ± 2.4)% of the 
J/ip counts in our sample come from ip' decays. (41.4 
± 0.9)% of these decays contain a neutral meson that 
decays into photons |45|. namely ip' ^ J/4' + tt^tt^, 
Ip' — >■ J/ip + 7r° and ip' J 14 + V- We will refer 
to these decay channels collectively as -0' J/ip -f 727. 
Simulations show that most of the decays into neutral 
mesons are either not detected in the central arm accep- 
tance or are rejected by the 7 energy cut, leaving an esti- 
mated 6-20 counts in the low mass distribution of e^e~ 7 
(Fig. [TT|) . Contributions from ip' — > 7 + Xc — > 27 -t- J/-0 
decays are expected to be no larger than three counts. 

The contribution from B decays in the e^e~ 7 sam- 
ple was calculated using the bottom cross section mea- 
sured by PHENIX [36] . The contribution of B decays to 
J/ 4 plus at least one photon is less than 3 counts in the 
entire e^e~7 sample. 

The number of Xc decays was obtained by fitting the 
background and the simulated Xc line shapes to the mea- 
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FIG. 10: (Color online) (a) Simulated e^e~ -I- 7 invari- 
ant mass distribution from pythia events containing pri- 
mary J/xp decays. The line is the combinatorial back- 
ground obtained using simulated mixed events from the 
same sample (top), (b) The e"'"e~7 mass distribution in 
the data, after combinatorial background subtraction, where 
Mg+g- [GeV/c'^ ] <2.9. The lines are empirical fits as ex- 
plained in the text. Note that the simulated ip' J/ip -|- 7 
shape is arbitrarily normalized. 
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FIG. 11: (Color online) Mass (a) and 7 energy (b) distri- 
butions of ip' J/ip -|-neutrals — > e^e~ 7 obtained from 
simulations. Dashed lines in panel (b) represents the photon 
energy cut applied in this analysis. Appropriate scaling indi- 
cates that such events contribute between 6 and 20 counts to 
the correlated e^e~ 7 distribution. 



sured e+e~7 mass distribution (Fig. [T2t . The back- 
ground includes two sources: the mixed event back- 
ground from random e'^e" + 7 combinations, and the 
correlated background discussed previously. The corre- 
lated background was fit to a Gaussian and a Landau 
distribution where the maximum of the correlated back- 
ground was set by the photon energy cut. In addition, 
the '0' ~^ J /4 + '^7 background was used as a third 
shape in estimating the systematic error. However, it 
must be emphasized that this background cannot explain 
the magnitude of the correlated background. The vari- 
ations introduced by using the three distributions con- 
tribute to the type B systematic errors. The fitting 
parameters included the combinatorial background nor- 
malization, the amplitude of the correlated background 
and, when used, the width of the Gaussian and Lan- 
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FIG. 12: (Color online) Top plot: e'^e~j invariant mass dis- 
tribution where the dielectron is required to have an invariant 
mass within J/tp mass region. The like-sign dielectron distri- 
bution is subtracted as described in the text. Bands represent 
the contributions from various sources: solid line - Xc signal; 
long dashed line - mixed event uncorrelated background; short 
dashed line - correlated background as described in the text. 
The correlated background was fit to a Gaussian for this plot. 
The widths of the bands reflects the Icr variations in the fit. 
The bottom plot is the Xc signal after subtraction of the back- 
grounds. 



dau shapes (the ip' — J/ip + nj shape was fixed from 
simulations), and the normalization of the simulated 
Xc mass distribution. The fitted Xc mass spectrum re- 
turned an average value 96 ± 24 counts in the Xc mass 
range Afg+e-^ — M^+^- e [0.3,0.6] GeV/c^ when fitting 
the three different line shapes to the correlated back- 
ground. The signal/background, including the correlated 
background, was 1/5. The number of Xc counts changed 
by ± 4.6% when using different line shapes for the corre- 
lated background (Gaussian, Landau or ip' J/tp +nj 
shapes) . 



several px bins. While the statistical errors are large, 
the dependence of the acceptance on the pt of the Xc is 
mild, hence the error in the mean conditional acceptance 
is small. We obtain (e^Jej/^) = (12.0 ± 0.4) %. 

Tests of the fitting procedure and the conditional ac- 
ceptance calculation were performed using several dif- 
ferent simulated data sets with varying amounts of 
Xc signal, J/ip , and backgrounds. The feed-down ob- 
served after full analysis of the six sets of simulated events 
correctly returned the fraction of Xc events with no sig- 
nificant bias. Variations in the minimum Ej criteria 
changed the measured feed-down in the simulation by 
1.7%. This variation is taken into account in the uncer- 
tainties as a type B error introduced by the analysis pro- 
cedure. When the photon energy resolution is changed 
in a manner consistent with the measured tt*' — >■ 2j 
mass resolution, both the conditional acceptance and the 
Xc counts returned from the fits change, leading to a vari- 
ation of the feed-down fraction by 1.6%. The list of all 
systematic uncertainties is shown in Table HIl 



TABLE II: Summary of the type B systematic uncertainties 
in the Xc feed-down fraction measurement. The total gives 
the sum of all errors in quadrature. 



syst uncertainty 


contribution 


type 


7 ID 


0.7% 


B 


energy resolution 


1.6% 


B 


Xc polarization 


1.8% 


B 


correlated background line shape 


1.5% 


B 


J/t/i continuum 


0.1% 


B 


fit procedure 


1.7% 


B 


Xc momentum dependence 


1.1% 


B 


TOTAL 


3.6 % 





The final Xc feed-down fraction using ([T2|) is 



F^/'f' = 32 ± 9% 

AC 



(13) 



when taking the quadratic sum of the statistical and sys- 
tematic uncertainties. 



D. Feed-down fraction result 

The fraction of J/ip counts coming from Xc decays is 

To find the mean conditional acceptance, (^e^^/sj/^'), the 
conditional acceptance shown in Fig. [HI must be convo- 
luted with the Xc Pt distribution. An estimate of the 
Xc Pt distribution was obtained by fitting a two dimen- 
sional e^e~7 mass vs. pt distribution to a Xc signal plus 
backgrounds and extracting the number of Xc counts in 



V. J/tP ANALYSIS IN THE FORWARD 
RAPIDITY REGION 

This section describes the analysis performed to ob- 
tain the inclusive J/ijj dimuon yield at forward rapidity 
1.2 < \y\ < 2.4. Section |V3 describes the J /ip signal 
extraction from the dimuon spectrum and related uncer- 
tainties. The response of the muon arm spectrometers 
to dimuon decays from the J/ip is described in section 
IV Bl Finally, the px and rapidity dependence of the 
J/ ip differential cross section and a summary of system- 
atic uncertainties is reported in section IV CI 
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FIG. 13: (Color online) Invariant mass distribution of 
dimuons in the J/ip mass region. The components of the 
spectrum are the combinatorial background estimated from 
a mixed-event technique, an acceptance modified (Fig. I14|) 
two-Gaussian J/i/) signal, and an acceptance modified expo- 
nential continuum (Eq. llSp . 



sian which describes the hne shape of J/ tp in the Monte 
Carlo, and acceptance dependence: 



Acc {M^^) 



cont^ 



(15) 



{l-fG2)G{M^^,Mj/^,aGi) 
/g2G {Mfj,fj,,Mj/^ + 5m,(JG2) 

1 (M^^-M)^ 



27r(T 



where Acc{M^^) is the mass dependence of the dimuon 
acceptance in the rapidity 1.2< \y\ <2.4 estimated using 
dimuon simulation (Fig. [T4|) . A^p is the amplitude of the 
J/f/; signal with mass Mj/^ composed of a Gaussian of 
width (7(31 and a second Gaussian of width aQ2 shifted by 
5M in mass. /g2 describes the fractional strength of the 
second Gaussian. The normalization of the continuum 
contribution is Acont and its exponential slope is b~^^^. 



A. J /ip fJ-^ 1^ signal extraction 

The dimuon invariant mass spectrum was obtained 
from the muon sample selected according to the cri- 
teria described in Sec. HH The MuID trigger condi- 
tion is emulated offline. In order to make sure the real 
J/ip candidate fired the MuID trigger, at least one muon 
of the dimuon pair is required to match a road from the 
trigger emulator. 

The decomposition of the dimuon background is very 
similar to that described in Sec. IIII Al for dielectrons. 
The combinatorial background was estimated using the 
mass spectrum of random pairs formed by pairing op- 
posite sign muon candidates from different events. The 
muons of the mixed pair are required to have vertices 
that differ by no more than 3 cm in the beam direction. 
The mixed event spectrum was normalized by the factor 



( Tysame \ / Tysame 
^ /^+ f-L~ fj.~ 



I ( AjraixedX ( i\Tmixed\ 



(14) 



^ . are the number of pairs formed 



where iV,':™'= and NZ"""^ 

from two muons in the same or in mixed events, re- 
spectively. The mass spectrum of the dimuons in the 
J/ tp mass region is shown in Fig. 1131 

The components of the correlated dimuon spectrum 
are muons sharing the same cc or bb ancestor, dimuons 
from Drell-Yan and the J/ijj and ip' resonances. There 
is no clean mass discrimination between the J/tp and 
ip' mass peaks in the muon arm spectrometers. How- 
ever the Ip' contribution is expected to be negligible in 
the peak integral compared to other uncertainties. The 
correlated dimuon mass distribution can be represented 
by a function i^(Af^^) including an exponential shape ac- 
counting for the continuum distribution, a double Gaus- 
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FIG. 14: (Color online) Mass dependence of the dimuon geo- 
metric acceptance in the muon arm spectrometers. 

The correlated mass distribution function F [M^^) was 
fit to the measured unlike-sign dimuon mass distribu- 
tion for each px and rapidity range using the maximum 
likelihood method. The combinatorial background, ob- 
tained from the normalized mixed event distribution, was 
also introduced in the fit with a fixed amplitude. The 
mass resolution obtained in the entire J /ip sample was 
(JGi/Mji^ =4% (tJGi = 125 MeV). The fitting parame- 
ters which determine the line shape of the J/ip peak (ctgi , 
f G2i <^M and /g2) obtained from the entire unbinned sam- 
ple were fixed when performing fits for individual pT and 
rapidity bins. The J /tp mass, Mjj,^^ was allowed to vary 
by 10% of its nominal value (3.096 GeV/c^ ) in the fit- 
ting procedure, the J jip and continuum amplitudes were 
constrained to avoid unphysical negative values, and the 
exponential slope was allowed to vary by 20% from a 
the value found in a fit to the entire (unbinned) sam- 
ple. For the systematic uncertainty evaluation, ]q2 was 
changed by 25% up and down, the fit was performed 
in two mass ranges: 1.8 < Af^^[GeV/c^ ] <7.0 and 
2.2< il/^^[GeV/c^ ] <6.0 and the combinatorial back- 
ground normalization a was varied by ±2%. Fig. [T3] 
shows the fitted function and its components for the 
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(]) projection on MuTr station 1 

FIG. 15: (Color online) Simulated (shaded area) and real 
data (bars) single muons distributed in the (j> coordinate of 
the MuTr. 

dimuon unlike-sign distribution for one of the rapidity 
bins. Two methods for counting the J/^ s were consid- 
ered: 1) using the fitted amplitude directly, or 2) from 
direct counting of dimuon pairs in the mass region 2.6 
< M^+^- [GeV/c^ ] < 3.6 with subtraction of the combi- 
natorial background and exponential continuum under- 
neath the peak in that same region. The standard devi- 
ations of the central values of the fits and of the signal 
extraction method variations are taken as type A sig- 
nal extraction systematic uncertainties, since these vari- 
ations are largely driven by statistical variations. The 
total number of J/ip counts was 16, 612±147''''^*± 112">"^' 
in the south muon arm and 16, 669 ± 145''*''* ± 115''^''* in 
the north muon arm. 



B. Di-muon acceptance and efficiency studies 

The response of the muon arm spectrometers to 
dimuons from J/ip decays was studied using a tuned 
GEANxS-based simulation of the muon arms and an of- 
fiine MuID trigger emulator. The MuID panel-by-panel 
efficiency used in these simulations was estimated from 
reconstructed roads in real data, or in cases with low 
statistics, from a calculation based on the operational 
history record for each channel. The MuID efficiency 
had a variation of 2% throughout the Run leading to a 
systematic uncertainty of 4% for the J/i/; yield. 

The charge distribution in each part of the MuTr ob- 
served in real data and the dead channels and their varia- 
tion with time over the run were used to give an accurate 
description of the MuTr performance within the detector 
simulation. The azimuthal distribution of muon candi- 
dates in real data and simulated muons from J /tp decays 
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FIG. 16: Transverse momentum dependence of the north 
and south average muon arms acceptance x efficiency for 
J/t/i dimuon decays in 1.2 < \y\ < 2.4. Shaded bands are 
the uncertainties of the estimates described in the text. 

using the pythia simulation are shown in Fig. [15] The Z 
vertex distribution of simulated J/ip decays is the same 
as that observed in real data. The px distribution in 
the MuTr obtained in simulation was also weighted ac- 
cording to that observed in the real data. The relatively 
small differences in the real and simulated (p distributions 
(Fig. [T5]) are thought to be due primarily to missing 
records for short periods of time in the dead HV channel 
records. These differences are estimated to change the 
J/4> dimuon yields by up to 6.4(4.0)% in north(south) 
arms. Run- by-run variations of the MuTr single muon 
yields are estimated to affect the final J/ip yields by an 
additional 2%. 

The J/ip acceptance x efficiency {Ae) evaluation used 
a PYTHIA simulation with several parton distributions as 
input to account for the unknown true rapidity depen- 
dence of the J/ tp yield leading to variations of 4% in the 
final acceptance. Fig. [Tni shows the overall pr depen- 
dence of Ae for J/ip dimuon decays. The uncertainties 
related to the knowledge of the detector performance are 
point-to-point correlated between different pT and differ- 
ent rapidity bins. The uncertainty in the dimuon accep- 
tance caused by lack of knowledge of the J/ tp polarization 
was studied using the detector simulation. The first re- 
sults in PHENIX at forward rapidity [1^ indicate that 
the J/V' polarization is no larger than ±0.5 for px <5 
GeV/c (in the Helicity frame). For this polarization vari- 
ation, the simulations show one standard deviation vari- 
ations between 2% and 11%, with the largest variation 
occurring for px <l GeV/c and y ~ 1.2. For px > 5 
GeV/c , where there are no polarization measurements 
we consider polarizations anywhere between ±1, and find 
variations no larger than 5%. These deviations are con- 
sidered as type B uncertainties. 



C. J/tp dimuon cross section result 

The differential cross section for each px bin was calcu- 
lated according to Eq. The systematic uncertainties 
involved in this calculation are listed in Table UlIl 
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TABLE III: List of the systematic uncertainties in the 
J/t/i dimuon yield measurement. Ranges indicate pt depen- 
dence. 



description 


relative uncertainty 


type 


signal extraction 


1.8% - 35% 


A 


MuID efficiency 


4% 


B 


MuTr acceptance 


6.4%(north), 4.0%(south) 


B 


run-by-run fluctuation 


2% 


B 


Monte Carlo J /ip input 


4% 


B 


J/'4> polarization 


2% - 11% 


B 


luminosity 


10% 


C 
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FIG. 17: Transverse momentum dependence of the 
J/xp dimuon differential cross section obtained in the muon 
arms in 2006 and 2008 Runs. 



The differential cross section was independently ob- 
tained in the north and south muon arm spectrometers 
and for the 2006 and 2008 Runs. The measurements 
agree in all data sets for all pT points within one sigma 
statistical and systematic uncertainties. The averaging 
of these four momentum spectra is done using a weight 
for each data set based on the uncertainties for each that 
are uncorrelated between data sets. By definition the 
statistical and type A uncertainties are uncorrelated and 
while the type C is correlated. The uncertainties in the 
MuTr efficiency and run-by-run variations are also un- 
correlated between data sets. The MuID efficiency and 
the simulation input uncertainty are correlated between 
different spectrometer arms and run periods. Fig. [iTl 
shows the resulting average differential cross section for 
dimuons from J /ij} . The numbers are listed in Table [ 

The rapidity distribution was calculated as: 



N 



j CdtAe 



(16) 



where the number of J/i/" counts Nj/^ and the accep- 
tance X efficiency estimates were performed for each ra- 
pidity bin. All these results are shown in Fig. [18] and the 
numerical results are listed in Table IIXI 



VI. RESULTS DISCUSSION 



This section presents a summary of the results re- 
ported in the previous sections and compares them with 
results obtained in other experiments as well as predic- 
tions from several different production mechanism calcu- 
lations. The rapidity dependence of the J /ip yield is com- 
pared to models using various parton distribution func- 
tions (PDF) in Sec. IVI Al The total J/-0 cross section 
is derived from the rapidity distribution and discussed in 
Sec. IVIBI The J/i/' differential cross section dependence 
on Pt is compared to empirical scaling laws observed at 
lower energies as well as different charmonium hadroniza- 
tion models in Sec. IVI CI The measured fraction of the 
J I ip yield coming from ip' and Xc decays is compared to 
other experiments in Sec. IVIDI The consequences of the 
results presented in this article on recent charmonium 
measurements in p{d)+A and A-l-A collisions is the sub- 
ject of the Sec. IWEl 

The models used in our comparisons were described in 
Sec. m namely, the Color Evaporation Model (CEM), the 
Color Singlet Model (CSM) and Non-relativistic QCD 
(NRQCD). The CEM used FONLL calculations for the 
charm cross section and CTEQ6M as the parton dis- 
tribution function pol. Iso} . For the CSM comparison, 
we used the recent NLO calculation only for the direct 
J/i{j yield at RHIC energy and PHENIX rapidity cov- 
erage We used two NRQCD calculations in our 
comparisons. The calculation performed for the direct 
J/ ip plus Xc feed-down in [20|] uses NLO diagrams for the 
color singlet and color octet states with a long range ma- 
trix element tuned from experimental hadroproduction 
[2l| and photoproduction 22, 23] results. This calcula- 
tion is only available for the differential px dependent 
cross section. An older calculation, performed for the 
same direct J/tp plus Xc feed-down with LO diagrams 
5l[ , also provides the rapidity dependence and total cross 



sections for different PDFs. No similar attempt has been 
made with the new calculations. The differential pt de- 
pendent cross section calculation involves the emission 
of a hard gluon which determines the shape of the char- 
monium Pt spectrum. The amplitude of the hard gluon 
emission cannot be calculated for pT < 2 GeV/c because 
of infrared divergences. This problem is circumvented in 
the older calculation by empirically constraining the low 
Pt nonperturbative soft gluon emission to obtain the ra- 
pidity dependence, dcr/dy. In both NRQCD calculations 
there is a prevalence of color octet states in the direct 
J 1^ contribution. 
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FIG. 18: (Color online)Rapidity dependence of the J/ij) yield 
combining dielectron < 0.35 - full squares) with dimuon 
channels (1.2 < \y\ < 2.4 - full circles) along with the fits used 
to estimate the total cross section. Lines correspond to the 
three fitting functions described in the text. Also shown are 
arbitrarily normalized model predictions (NRQCD [5j| , CEM 
[4a,i5Q] and CSM jj]). 



A. J/ip Rapidity dependence 

The rapidity distribution of the J/^ dilepton cross sec- 
tion is shown in Fig. [TSland in Table HXl The data points 
are grouped into three rapidity ranges, corresponding to 
the different detectors used in the measurement: south 
muon arm (—2.4 < y < —1.2), central arms (|?;| < 0.35) 
and north muon arm (1.2 < y < 2.4). The systematic 
uncertainties represented by the boxes are point-to-point 
correlated for data points in the same group and are un- 
correlated between different groups. All points have a 
global uncertainty of 10% coming from the minimum- 
bias trigger efficiency estimate. 

In order to compare the shape of the rapidity distribu- 
tion, we normalized the CEM, CSM and NRQCD predic- 
tions to the integral of the measured data in FigllHl All 
models use the CTEQ6M PDFs. The NRQCD model is 
also available with the GRV98 and the MRST99 PDFs. 
The theoretical rapidity distributions exhibit a similar 
shape when using CTEQ6M. A very different rapidity 
distribution is obtained when the NRQCD prediction is 
calculated using GRV98 and MRST99 (MRST99 is not 
shown in the figure). These observations suggest that 
the choice of PDF plays the most important role in de- 
scribing the shape of the J/ip rapidity distribution. The 
rapidity shape also appears to be independent of the feed- 
down contributions, since the CSM has a similar shape 
to the CEM and NRQCD model, despite the fact that it 
contains only direct J/^ contributions. The PDF which 
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FIG. 19: (Color online) Energy dependence of the parameter 
c where da/dxp = (1 — xfY is fitted to xf distributions 
of J/ip production in fixed target experiments [s^ - lsol and 
in the PHENIX rapidity distribution. The parameters (a = 
27.8 ± 1.5, 6 = (143 ± 11) GeV) are obtained from a fit to the 
experimental data points. 



best describes the data is CTEQ6M and we use this for 
the remaining comparisons. 

An empirical description of the J /ip yield used in some 
fixed-target experiments with large coverage is based on 
the Feynman xp form [s^ . 



da 



A{i-\xp\r. 



(17) 



We can convert to a rapidity distribution by writing 
da dxp 



da 

dy dxp dy 
« 2A{1-\xf\ 



(18) 



'(pi 



coshy 



where px — 1.73 GeV/c is the average of the J/ip px dis- 
tributions over all measured rapidities. The fit returned 
c — 16.3 ± 0.4 with probability of 31%, where statisti- 
cal and systematic uncertainties are summed in quadra- 
ture(Fig. [T8| . Fig. [T9l shows that c scales approximately 
a.s c = a/ (1 + b/y/s). This extrapolation of the rapidity 
dependence can be used to estimate the total cross sec- 
tion from measurements with limited rapidity coverage 
and will be used as one method to calculate the total 
cross section from the present measurement. 



B. Total cross section of inclusive J/i/) 

The total cross section was estimated from different 
empirical functions fitted to the rapidity distribution - a 
double Gaussian, the xp scaling function ^TE\\ described 
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above, and a "Fermi" function: 



Rapidity distributions based on charmonium production 
models were not used in the total cross section in order 
to avoid any theoretical bias. 

The correlated uncertainties between data points mea- 
sured in each spectrometer were propagated to the fit 
uncertainty by allowing the points to move coherently in 
the rapidity range covered by that spectrometer. Table 
II VI shows the total dilepton cross section and the prob- 
ability for each function used in the fit. The final cross 
section is obtained from the average of the numbers from 
each fit function weighted according to their ^ prob- 
ability. The systematic uncertainty from the unknown 
rapidity shape is taken from the standard deviation be- 
tween the three fitting functions. Based on these fits, we 
conclude that the PHENIX rapidity acceptance covers 
56 ± 2 % of the total J/?/' cross section. The J/'0 cross 
section reported in this paper is 180.7±2.0''''^*±ll^y'''n6, 
in agreement with our previous result with a reduction 
in the statistical and systematic uncertainties. 



TABLE IV: Estimate of the total dilepton J/i/) cross section 
from the three fitting functions, together with the weighted 
average and a comparison to the result obtained in our 
previous measurement. The measured total cross sections 
have an additional 10% global uncertainty. 



estimating function 


prob. 




Biicyji^ (nb) 


xf scale fen, EqIISI 


0.30 


170.8 




double Gaussian 


0.79 


183.5 


_|_ gstat _|_ -j^-j^syst 


Fermi fen, EafT9l 


0.70 


182.0 




AVERAGE 




180.7 


± 2.0"''" ± 12"^"' 


2005 Run result [10] 




178 ± 


3 gstat _j_ ggsyst 



Table IVl presents the measured total J/ip cross section 
and the expectations from the three production mod- 
els considered in this text. The experimental direct 
J/ip cross section is estimated assuming that the feed- 
down fraction of Xc and ip' measured at midrapidity is 
the same at forward rapidity. The feed-down from B 
mesons is only significant at high pT and is not consid- 
ered in the total cross section. The total cross section es- 
timated using the CEM is the only one which agrees with 
the experimental result, although the cross section calcu- 
lation includes the scale factor (Sec. H]) obtained from 
J/ip measurements. The NRQCD includes color singlet 
and color octet states, and as mentioned at the begin- 
ning of this section, cannot be extrapolated to low px to 
obtain the rapidity distribution without the addition of 
an empirical constraint. 



TABLE V: Comparison of the measured J/t/j cross section 
with the three models considered in this text. Direct 
J/t/i cross sections are obtained assuming that the Xc and 
ip' feed-down fractions measured at midrapidity are the same 
at forward rapidity. Type A, type B and type C errors are 
quadratically summed in the measured result. 





direct J/t/j 


inclusive 


CEM 




169 ± 30 nb 


NLO GSM 


53 ± 26 nb 




LO NRQGD 




140 ± 5 nb 


Measured 


105 ± 26 nb 


181 ± 22 nb 



C. J/tp pt distribution 

The Pt -dependent dielectron differential cross section 
at midrapidity is compared to other p+p and p + p ex- 
periments in Fig. BUT aV The shapes of the transverse 
momentum distributions follow the well known "ther- 
mal" exponential behavior for pr < 2 GeV/c and a hard- 
scattering power law behavior at high px ■ The hard pro- 
cess scales with xt — 2pt/\/s {^/s^ E(fia / (Pp — G{xt)) 
(60| for all collision energies, as can be seen in Fig. EOfb). 
where n — 5.6 ± 0.2 n is related to the number of 

partons involved in the interaction. A pure LO process 
leads to n = 4, hence, NLO terms may be important in 
J/ip production. |62H65j 

The Pt dependence of the J /ip differential cross sec- 
tions measured at forward and midrapidity are shown 
in Fig. [21] along with theoretical calculations where the 
absolute normalization is determined in the calculations. 
The CEM and the NRQCD (for pr > 2 GeV/c ) provide 
reasonable descriptions of the pr distribution, whereas 
the CSM disagrees in both the normalization and the 
slope of the px distribution, indicating that NLO color 
singlet intermediate states cannot account for the direct 
J /ip production. However, the NLO CSM calculation 
gives a good descrigtion of the J /ip polarization mea- 
sured by PHENIX [l3,S|. Attempts are being made to 
extend the CSM to NNLO. Preliminary NNLO CSM cal- 
culations performed for pT > 5 GeV/c [13 shows a large 
increase in the yield, but still under- predict the experi- 
mental results. None of these theoretical models consider 
the B-meson decay contribution to the J/ip yield. The 
fixed-order plus next-to-leading-log (FONLL) fs'] calcu- 
lation of these decays is also plotted in Fig. [51] and 
has a reasonable agreement with STAR measurements 
using J/ip -hadron correlations f6TI|. According to this 
calculation, the i?-meson contribution to the measured 
J/ip inclusive yield is between 2%(1%) at 1 GeV/c and 
20%(15%) at 7.5 GeV/c in the mid(forward)-rapidity re- 
gion with large theoretical uncertainties. 

The Pt dependence of the J/ ip yield is harder at midra- 
pidity, as seen from the ratio between the forward and 
midrapidity differential cross sections versus pt shown 
in Fig. I^TTbottom). The figure also includes the for- 
ward/midrapidity yield ratios from the theoretical mod- 
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FIG. 20: (Color online) pr distribution (a) and xt distribution (b) of J/iP yield in PHENIX, STAR |6i], ISR UAl ^ 
and CDF[i,[2l| at y ~ 0. 



els using their mean values and assuming that theoretical 
uncertainties in these ratios cancel out. All of the mod- 
els predict a downward trend, but the CEM and NRQCD 
calculations do not follow a slope as large as the data. 

The mean transverse momentum squared {px} was cal- 
culated numerically from the pt distribution. The cor- 
related uncertainty was propagated to {pt} by moving 
\ow-pT and high-p^ data points coherently in opposite 
directions according to their type B uncertainty. The re- 
sults with the propagated type A and type B uncertain- 
ties are listed in Table IVll The table also contains (px) 
for Pt < 5 GeV/c for a direct comparison with previous 
PHENIX results 68]. As expected, the mean transverse 
momentum squared at midrapidity is larger than at for- 
ward rapidity. 



TABLE VI: Mean transverse momentum squared in 
(GeV/c )^ of J/ip and ip' for different rapidity andpT ranges. 
Uncertainties are type A and type B respectively. 



J/ij 1.2 < \y\ < 2.4 3.65±0.03±0.09 3.45±0.03±0.08 
J/V) \y\ < 0.35 4.41±0.14±0.11 3.89±0.11±0.Q9 

^' \y\ < 0.35 4.7 tli, ±0.2 A.7tli, ±0-2 



D. Charmonia ratios and J/tp feed-down fractions 

The transverse momentum dependence of the 
4>' /{J/'ip ) yield ratio (Fig. [51 bottom) is consistent 
with that observed in other experiments. Fig. [22] shows 
the collision energy dependence of the i/;' /(J/V' ) yield 
ratio in light fixed target experiments and p+p or p + p 
colliders. In this figure, the ratios from p+p experiments 
were calculated using the reported J/^p and ip' cross 
sections for px > 5 GeV/c together with their point-to- 
point uncorrelated uncertainties'^. The B meson decay 
contribution was removed from the J/'ijj and tjj' yields, in 
the case of the CDF experiment. Only E705 has broad 
coverage (—0.1 < a;^ < 0.5). The other experiments in 
this figure have a rapidity coverage of \y\ < 0.6. A weak 
trend of increasing ip' /{J/tp ) yield ratio for higher 
collision energy (Fig. [51]) and for higher pT (Fig. [S]) can 
be observed. As mentioned earlier, the ip' feed-down 
fraction of (9.7 ± 2.4)% is in agreement with the world 
average of (8.1 ± 0.3)% calculated in 69]. 

The feed-down fraction obtained from our 
Xc ^ J/ip + 7 measurement is compared with other 
experiments over a broad range of collision energy and 
xp, as well as over many different colliding species (Fig. 

[23| . The value measured in this work, F^/''' = (32±9)%, 



^ This may be an overestimate of the systematic errors, given that 
a good fraction of the J/tp and ip' yields may be correlated. 



21 



p+p^ J/\|/+X \fs=200GeV 

- - NLO CSM (direct J/\|/) 

— NLONRQCD 
-CEM 

--- B-meson— > J/\j/+X 
. lykO.35 X 1000 
1.2<lyl<2.4 




[GeV/c] 



FIG. 21: (Color online) Top: Transverse momentum depen- 
dence of J/i/) yield in \y\ < 0.35 and 1.2 < \y\ < 2.4 along with 
predictions based on CSM [l3|, NRQCD CEMIS, H^l 
and B-meson decay based on FONLL calculation All 
models use CTEQ6M. Theoretical uncertainties are repre- 
sented as bands. Note that the midrapidity points are scaled 
up by a factor of 1000. Bottom: Ratio of the forward and 
central rapidity pr spectra and corresponding theoretical pre- 
dictions. 
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FIG. 22: (Color online) Collision energy dependence of the 
■i/)' /{J /tp) dilepton cross section ratio obtained in p-\-p and 
light fixed target p-l-A collisions [3, 0, [S^, ItoI - ItI ] compared 
to the CEM estimate [4^. Statistical and systematic errors 
were quadratically summed. 
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is consistent with the world average of (25 ±5)% after 
accounting for A dependencies in the fixed target 
experiments j69|.^ 

Combining the results of feed down from the ip' and the 
Xc we obtain a total J feed-down fraction measured in 
the midrapidity region of (42 ±9)%. 

E. Outcomes for heavy ion collisions. 

The feed-down fractions from the and Xchave im- 
portant implications for survival rates of charmonium 
states when either ordinary nuclear matter or high en- 
ergy density nuclear matter is involved. Because of their 
larger size compared to the J/ij} , excited charmonium 
states may have a different breakup cross section in nu- 
clear matter. This effect can modify the feed-down frac- 
tions in p+A collisions. On the other hand, if the cc is 



* The world average was obtained in after extrapolating the de- 
pendence of the estimated path length in nuclear matter for the 
p+A fixed target experiments. 



FIG. 23: (Color online) Collision energy dependence of 
Yc feed-down to J/'i/> measured in p+p and p-l-A collisions 
[toI. It^-IsoI] compared to the CEM calculation [ii]. 

not formed as a color singlet, it can cross the nuclear mat- 
ter as a colored preresonant state [Slj. If this were true, 
the breakup cross section of J/^ and t/'' should be the 
same and there would be no modification of the tp' feed- 
down fraction in p-l-A collisions, whereas a possible mod- 
ification can occur for the Xc since it is expected to be 
formed mainly as a color singlet. Given the large sta- 
tistical uncertainties in all measured tp' / {J /ip ) ratios 
shown in Figs. H] and [521 differences between p + p and 
p+K are impossible to see, and therefore no conclusion 
about possible cold nuclear matter effects can be made at 
this time. The same is true for the Xc feed-down fraction 
in Fig. 1231 Higher precision measurements of charmo- 
nium states in p+p and d-l-Au collisions in the future 
may allow an improved determination of these possible 
cold nuclear matter effects on the feed-down fraction. 

The behavior of charmonium states in the high den- 
sity, hot nuclear matter created in heavy ion collisions, 
has long been of interest [l'] . Spectral function computa- 



22 



tions [3] indicate that the Xc and if)' states should disso- 
ciate at a lower temperature in hot nuclear matter, due 
to color screening, than the J/V'. One of the most im- 
portant implications of the observed feed-down fractions 
is that the complete dissociation of the Xc and states 
would lead to a (42 ± 9)% J/tj: suppression. The mea- 
sured nuclear modification factor of J /ip mesons in cen- 
tral Au-I-Au collisions at -^5 — 200 GeV [l^l implies a 
suppression of (74 ± 6)% at midrapidity and (84 ± 6)% 
at forward rapidity. Hence, the complete dissociation of 
the excited states of charmonium and the associated loss 
of the J /ip yield cannot completely explain its suppres- 
sion observed in Au-|-Au collisions. Cold nuclear mat- 
ter effects and the possible dissociation of direct J/ ^ by 
color-screening could presumably account for the remain- 
ing suppression. 



VII. SUMMARY AND CONCLUSIONS 

In conclusion, we have measured the yields of the three 
most important charmonium states in p+p collisions at 
— 200 GeV, where gluon fusion is expected to be the 
dominant production process. The rapidity dependence 
of J/ip supports the use of CTEQ6M to describe the 
gluon distribution in protons. The inclusive J /ip yield is 
in agreement with current models which involve a initial 
formation of colored charmonium states, as in the GEM 
or the color octet states of the NRQGD models. The in- 
clusive J/ ij} yield observed at midrapidity is composed of 
9.6 ± 2.4% of V' decays and 32 ± 9% of Xc decays. This 
result is in agreement with what was observed in other 
experiments. Given the current large statistical uncer- 
tainties, no conclusion can be made about collision en- 
ergy or pt dependence of these fractions. Finally, this 
J cross section measurement and feed-down fractions 
will play an important role in current studies of cold nu- 
clear matter and the hot, dense matter formed in heavy 
ion collisions. 



Acknowledgments 

We thank the staff of the Gollider-Accelerator and 
Physics Departments at Brookhaven National Labora- 
tory and the staff of the other PHENIX participating 
institutions for their vital contributions. We also thank 
Jean Philippe Lansberg, Mathias Butenschon and Ra- 
mona Vogt for the valuable GSM, NRQGD, GEM and 
FONLL calculations in the PHENIX acceptance. We ac- 
knowledge support from the Office of Nuclear Physics in 
the Office of Science of the Department of Energy, the 
National Science Foundation, a sponsored research grant 
from Renaissance Technologies LLG, Abilene Ghristian 
University Research Gouncil, Research Foundation of 
SUNY, and Dean of the GoUege of Arts and Sciences, 
Vanderbilt University (U.S. A), Ministry of Education, 
Gulture, Sports, Science, and Technology and the Japan 



Society for the Promotion of Science (Japan), Gonselho 
Nacional de Desenvolvimento Gientffico e Tecnologico 
and Fundagao de Amparo a Pesquisa do Estado de Sao 
Paulo (Brazil), Natural Science Foundation of Ghina 
(P. R. Ghina), Ministry of Education, Youth and Sports 
(Gzech Republic), Gentre National de la Recherche Sci- 
entifique, Gommissariat a I'Energie Atomique, and In- 
stitut National de Physique Nucleaire et de Physique 
des Particules (France), Ministry of Industry, Science 
and Tekhnologies, Bundesministerium fiir Bildung und 
Forschung, Deutscher Akademischer Austausch Dienst, 
and Alexander von Humboldt Stiftung (Germany), Hun- 
garian National Science Fund, OTKA (Hungary), De- 
partment of Atomic Energy and Department of Sci- 
ence and Technology (India), Israel Science Foundation 
(Israel), National Research Foundation and WGU pro- 
gram of the Ministry Education Science and Technol- 
ogy (Korea) , Ministry of Education and Science, Russian 
Academy of Sciences, Federal Agency of Atomic Energy 
(Russia), VR and the Wallenberg Foundation (Sweden), 
the U.S. Givilian Research and Development Foundation 
for the Independent States of the Former Soviet Union, 
the US-Hungarian Fulbright Foundation for Educational 
Exchange, and the US-Israel Binational Science Founda- 
tion. 



Appendix: Data Tables 



TABLE VII: J/ip differential cross section in the midrapid- 
ity region (|y| < 0.35) followed by point-to-point uncorre- 
lated (statistical and uncorrelated systematic uncertainties) 
and correlated systematic uncertainties. Global uncertainty 
is 10%. 



PT [GoV/c ] 


1 




[n&/(GeV/c f] 






dydprp 






value 


uncor. 


corr. 




0-0.25 


4.9 


0.5 


0.6 




0.25-0.5 


3.9 


0.3 


0.5 




0.5-0.75 


3.9 


0.2 


0.5 




0.75-1 


3.5 


0.2 


0.4 




1-1.25 


3.19 


0.17 


0.38 




1.25-1.5 


2.21 


0.14 


0.27 




1.5-1.75 


1.69 


0.12 


0.2 




1.75-2 


1.42 


0.1 


0.17 




2-2.25 


95 


8 


12 


XlO"^ 


2.25-2.5 


66 


7 


8 


XlO"^ 


2.5-2.75 


56 


6 


7 


XlO"^ 


2.75-3 


37 


5 


5 


xlO"^ 


3-3.25 


29 


4 


4 


xlO"^ 


3.25-3.5 


19.9 


+ 3.5 
-3.6 


2.4 


xlO^^ 


3.5-3.75 


13.6 


2.8 


1.7 


xlO"^ 


3.75-4 


10.6 


2.5 


1.3 


xlO"^ 


4-4.25 


6.89 


+ 2.3 
-2.2 


0.8 


xlO"^ 


4.25-4.5 


4.4 


+ 1.7 
-1.6 


0.5 


xlO"^ 


4.5-4.75 


4.7 


+ 1.5 
-1.4 


0.6 


xlO"^ 


4.75-5 


3.1 


1.2 


0.4 


xlO"^ 


5-6 


1.35 


+ 0.31 
-0.32 


+0.17 
-O.IS 


xlO"^ 


6-7 


4.1 


1.3 


+ 0.5 
-0.6 


xlO"^ 


7-8 


1.6 


0.7 


0.2 


xlO~^ 


8-9 


0.37 


+ 0.37 
-0.22 


+ 0.05 
-0.06 


xlO"^ 
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TABLE VIII: V' differential cross section at |77| <0.35 fol- 
lowed by point-to-point uncorrelatod (statistical and uncor- 
rclatcd systematic uncertainties) and correlated systematic 
uncertainties. (41()l)al uueertaiiitv is !()'/(. 



PT [GeV/c ] 



[pb/(GeV/c f] 



value 



0- 1 

1- 2 

2- 3 

3- 5 
5-7 
0-5 



67 
40 
15 
2.7 

95 



20 
11 
6 

+ 2-5 

<2.25 (90% CL) 
20 



+ 7 

-6 

3 



+ 17 
-15 



TABLE IX: Rapidity dependence of J/'ip yield followed by 
point-to-point uncorrelated (statistical and uncorrclated sys- 
tcniatic uncertainties) and correlated systcnnatic uncertain- 
ties. (Tli)bai uueertaiiitv is lO'X. 



rapidity 




Si! <i« 






value 


uncor. 


corr. 


-2.325 


10.9 


1.9 


1.0 


-2.075 


17.6 


0.5 


1.5 


-1.825 


24.4 


0.4 


1.9 


-1.575 


31.5 


0.5 


2.2 


-1.325 


41.2 


1.1 


5.3 


-0.3 


49.0 


2.1 


5.4 


0.0 


45.6 


1.6 


5.0 


0.3 


46.1 


1.9 


5.1 


1.325 


40.7 


1.2 


5.7 


1.575 


33.6 


0.7 


3.0 


1.825 


25.6 


0.4 


2.5 


2.075 


18.9 


0.4 


1.9 


2.325 


13.9 


0.9 


1.4 



TABLE XL Di-muon J/ip yield in the forward rapidity region 
(1.2< \y\ <2.4) followed by point-to-point uncorrelated (sta- 
tistical and uncorrelated systematic uncertainties) and corre- 
lated systematic uncertainties. Global uncertainty is 10%. 



PT [GeV/c ] 






[n6/(GeV/c ) ] 












value 


uncor. 


corr. 




0.125 


3.49 


0.14 


0.50 




0.375 


3.28 


0.08 


0.49 




0.625 


2.85 


0.06 


0.40 




0.875 


2.43 


0.05 


0.16 




1 125 


2 04 


04 


18 




1.375 


1.57 


0.03 


0.16 




1.625 


1.194 


0.024 


0.12 




1.875 


85.1 


1.9 


7.9 


Xl0~^ 


2.125 


63.8 


1.6 


5.6 


X 10^^ 


2.375 


46.8 


1.2 


3.7 


XlO"^ 


2.625 


31.5 


0.9 


2.5 


xlO~^ 


2.875 


22.1 


0.7 


1.7 


xlO~^ 


3.125 


15.3 


0.6 


1.1 


xlO~^ 


3.375 


11.1 


0.5 


0.8 


xlO~^ 


3.625 


7.7 


0.4 


0.6 


xlO~^ 


3.875 


5.53 


0.27 


0.37 


xlO"^ 


4.125 


3.28 


0.21 


0.23 


xlO"^ 


4.375 


2.26 


0.16 


0.15 


XlO"^ 


4.625 


1.45 


0.13 


0.10 


xlO"^ 


4.875 


1.06 


0.11 


0.08 


xlO-^ 


5.125 


1.02 


0.10 


0.07 


xlO-^ 


5.375 


4.3 


0.6 


0.4 


xlO-^ 


5.625 


2.9 


0.4 


0.2 


xlO-^ 


5.875 


2.4 


0.4 


0.2 


xlO-^ 


6.25 


1.02 


0.15 


0.09 


xlO-^ 


6.75 


0.73 


0.13 


0.06 


xlO-^ 


7.5 


0.13 


0.04 


0.012 


xlO-^ 



TABLE X: V' / (•^/V' ) dielectron yield ratio measured at 
|r?| <0.35 followed by point-to-point uncorrelated (statisti- 
cal and uncorrelated systematic uncertainties) and correlated 

systematic uucertaiuties. 



PT [GeV/c ] 










value 


uncor. 


corr. 


0- 1 

1- 2 

2- 3 

3- 5 
5-7 
0-5 


1.69 
1.96 
2.3 
3.4 

2.1 


0.51 
0.53 
1.0 

+ 2.0 

<38 (90% CL) 
0.5 


+u.iy 

-0.11 
+ 0.23 
-0.15 
+0.5 
-0.3 
+0.5 
-0.4 
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